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ABSTRACT 


Two tracer dye tests were carried out on the North Saskatchewan 
River to define the hydraulic and mixing characteristics of the river 
reach between Edmonton and the Saskatchewan border for mid-to-late 
winter ice conditions. Discharges during the tests ranged from 108 m/s 
to al 53 m/s. Ice cover thicknesses varied between 0.3 m and 0.8 m 
although frazil accumulations approaching 2 m were found in some 
sections. The mean water velocities calculated from the travel time of 
the peak dye concentration was 0.65 m/s in the open water reach between 
Edmonton and Fort Saskatchewan, 0.55 m/s in the ice covered reach 
between Fort. Saskatchewan and Shandro and 0.49 m/s in the ice covered 
reach between Shandro and the Saskatchewan border. Composite hydraulic 
roughnesses varied between 0.016 to 0.038. These composite roughnesses 
can be used to extrapolate velocities at other discharges provided that 


the ice characteristics do not change significantly. 


Both linear and Fickian dispersion occurred in the test reaches. 
The linear dispersion parameter was found to be 0.00612 for the Edmonton 
injection and 0.00429 for the Shandro injection. These values are lower 
than the values measured on other ice covered rivers but this is due to 
the low values of hydraulic roughness found in these river reaches. The 
Fickian dispersion coefficient was found to be 344 m/s for the upper 
test reach and 291 m/s for the lower test reach. These values are 
consistent with values reported in the literature for channels of 
Similar width and roughness even with the presence of the ice cover. 
The location of the transition between linear and Fickian dispersion was 


found to depend on the position and number of injection points. 


Both the linear and Fickian dispersion models predicted much more 
symmetrical distributions than were actually observed. The models 
accurately predict dispersion in uniform channels but the variations in 
geometry in natural rivers causes additional mixing which the models do 
not predict. 
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1.0 INTRODUCTION 


1.1 Background 


Increased awareness of environmental issues has prompted a 
reevaluation of the impacts of effluent discharges on the water quality 
of many Alberta rivers. Water quality modeling techniques are now being 
used to reassess the impacts of effluents from existing sources as wel] 
as to assess the potential impacts of effluents from new sources. These 
models provide a better understanding of the effects of effluents on 
river water quality; however, the models can only be used confidently if 


the hydraulic and mixing characteristics of the rivers are known. 


The most critical conditions for water quality in rivers in 
northern Alberta usually occur during the winter period. In most of 
these rivers the lowest flows occur during the winter. These low flows 
reduce the effluent assimilation capacity of the rivers. As well, the 
ice covers which exist on these rivers during the winter block the 
reabsorption of oxygen, thus intensifying the effects of oxygen 
consumption by some effluents. Unfortunately, there is a distinct lack 
of information on the hydraulic and mixing characteristics for these low 


flow ice-covered conditions. 


Effluent discharges from the Edmonton area, such as the sewage 
treatment plant, are being reassessed using the water quality models. 
One of the first steps in this assessment process is to determine the 
hydraulic and mixing characteristics of the river system downstream from 
the effluent outfall during the critical low flow winter period. To 
this end, the Alberta Research Council was requested by Alberta 
Environment to conduct tracer dye studies on the North Saskatchewan 
River between Edmonton and the Saskatchewan border in March of 1990 to 
measure travel times and to determine the longitudinal dispersion 
characteristics in this river reach. 
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This report addresses the objectives of the tracer dye studies 
conducted on the -North Saskatchewan River in March, 1990. These 


objectives are to: 


d. determine the travel times relative to the hydraulic 
characteristics to extrapolate travel times over a range of 


possible discharges encountered under an ice cover, and 


a): evaluate the longitudinal dispersion characteristics relative 
to the hydraulic characteristics. 
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2.0 FIELD INVESTIGATIONS 
2.1 Study Reaches 


The 336.3 km length of the North Saskatchewan River between 
Edmonton and the Saskatchewan border was divided into two reaches for 
which separate dye tests were carried out. The length of the upstream 
reach was 198.1 km while the length of the downstream reach was 
189.6 km. The reaches overlapped for 51.4 km between Shandro and 
Duvernay to compare the velocity measurements at the beginning and end 
of a test. These reaches were also chosen so that each dye test would 


be approximately the same duration. 


Figure la-b illustrates the entire length of the river between 
Edmonton and the Saskatchewan border and indicates the pertinent river 
distances (from the mouth), the two study reaches, locations where dye 
was injected, and the sample sites. Table 1 summarizes some of this 
information. The dye was always injected at the upstream end of each 
reach and concentrations were measured at six locations downstream. The 
subreaches varied in length from 20 km to 50 km depending on the 
availability of access. Dye injections proceeded from downstream to 
upstream to eliminate background interference from residual dye 
remaining from upstream tests. 


2.2 Methodology 
Cxca ingjecuion 


For each dye test, a predetermined volume of 20% Rhodamine WT dye 
was mixed with an equal volume of methyl] alcohol. The methyl] alcohol 
was used to adjust the specific gravity of the dye mixture to that of 
water. Multiple injection points in the cross-section were used at both 
injection sites to reduce the initial concentrations and reduce the 
distance needed to complete the transverse mixing process. At injection, 
the flow width was determined and the injection points spaced evenly 
across the channel. Different procedures were used to inject the dye 
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Table 1 Summary of field program. 


Reach Length Date Number of Duration of Dye Volume 
(km) 1990 Sample Sites Test (hours) (liters) 
Shandro 
- Meridian Lege Mar .06-Mar.12 6 US 67-2 
Edmonton 
- Duvernay 198.1 Mar .15-Mar .20 6 123 Be) 0 


into the flow because the river was ice-covered at one site but ice free 
“at the other site. 


At the injection site 100 m downstream of the Shandro bridge the 
river was ice covered so a number of 20 cm diameter holes were augered 
through the ice to define the width of flow. The two injection holes 
were selected to divide the channel width into equal segments. Once 
the injection points were chosen, another series of holes were augured 
upstream to serve as a water supply to flush the injection apparatus. 
The injection apparatus consisted of two 1.5 m lengths of 8 cm diameter 
ABS pipe. 


The injection apparatus was lowered into the holes such that the 
outlets were located in the center of the flow, midway between the bed 
and the bottom of the ice cover. If a longer length of pipe was 
required to achieve the desired injection configuration, a 1m long 
section could be added to the original length of pipe by a friction 
coupling. Once the injection apparatus was properly positioned, it was 
held in place against the ice by a tripod located on the ice surface. 
The pre-mixed dye was then simply poured into the pipe via a funnel, 
followed immediately by 20 to 40 L of clean river water to ensure al] 
the dye was removed from the ABS pipe and from the container holding the 
dye. 


The Edmonton injection was done from the Cloverdale pedestrian 
bridge just downstream of the Low Level bridge. Approximately one-third 
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of the channel, on the inside of the bend, was blocked by ice so the 
three injection points were selected to divide the open water portion of 
the channel into equal segments. The injection was done by using ropes 
to lower three buckets, each with one-third of the dye, to points just 
above the water surface. Another set of ropes attached to the bottoms 
of the buckets were then used to invert the buckets, pouring the dye 
into, the; river. The buckets were then immediately dropped into the 
river to wash away any remaining dye. This procedure produces a surface 
injection; however, the momentum of the dye causes it to mix well with 


the water rather than float on the surface. 
2.2.2 Sampling 


The sampling sites were established at locations where access was 
available. Following injection at the upstream end of the river reach 
being characterized, two crews of two persons followed the dye cloud 
downstream, each of the crews sampling at alternating sites. At each 
site, a number of 20 cm diameter holes were augured through the ice to 
establish the flow depth and the ice thickness. The water level was 
referenced to benchmarks so that water surface slope profiles could be 
estab! ished. 


The dye cloud was sampled over three holes across the channel. The 
centerline sampling location was established in the main part of the 
flow and two secondary sampling holes were established midway between 
the centerline sampling hole and both the left and right banks. 
Sufficient samples were taken at each of the sampling sites to 
completely characterize the dye cloud. Two or more samples were taken 
to establish background fluorescence levels and at least 20 samples were 
taken over the duration of the dye cloud. Depending on the duration of 
the dye cloud the interval between samples varied between 5 minutes and 
2 hours. The secondary sample sites were established to assess the 
extent of lateral mixing and to improve the estimates of the total 
volume of dye passing the sampling site. 


As mentioned earlier, the entire dye cloud was measured as 
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completely as possible. Special care was taken to define the time of 
first rise and the time of peak. Generally, the frequency of sampling 
was increased during those two periods. On the falling limb of the 
concentration curve, sampling was continued until it was judged that the 
fluorescence was reduced to 20% of the peak and that the tail of the 
curve could be confidently extrapolated. In some cases, however, the 
tails of the dye cloud were cut short in order to give sufficient time 


to go to the next sampling site. 


As soon as the samples were collected, they were transported to the 
mobile laboratory which each crew had at its disposal. The samples were 
maintained at temperatures of about 0 - 2°C (very nearly that of the 
river) and immediately run through the Turner Designs Model 10 fluoro- 
meter and the results recorded. Where the scales on the fluorometer 
overlapped, readings were taken on both scales to ensure consistency 


over the range of the two scales of the fluorometer. 
2.2.3 Calibration and Data Reduction 


Two fluorometers were used during the field study - the ARC Turner 
Designs Model 10 field unit and the U of A Turner Designs Model 10 rack 
mounted unit. Each was kept with the same crew for the duration of the 
entire campaign. Prior to the field work, a set of standards was pre- 
pared and both instruments were calibrated at 20°C for the set of stan- 
dards over a number of scales. After the field work, the instruments 
were recalibrated both at 20°C for a number of standards over a number 
of scales, and at O°C for two standards over the two most used ranges. 
There was essentially no difference between the pre- and post-field work 
calibrations. 


A simple least squares linear relationship was established for each 
scale and the fluorescent values recorded in the field were transformed 
into true fluorescent values, assuming the samples were collected at 
200% The final adjustment was made to compensate for temperature. 
Each value was multiplied by 0.62, the ratio of fluorescence between 
20°C dnd 7G. The data was then plotted and the pertinent character- 
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istics computed. 


2.3 Discharge Measurements 


Discharge measurements necessary to characterize the hydraulic 
conditions of each of the reaches were made by the Water Survey Section 
of Alberta Environment and the Water Survey of Canada. Three discharge 
measurements were generally made in each test reach. The first 
measurement was made as close as possible to the time and location of 
the injection, the second measurement was made at a sample site near the 
midpoint of the reach, and the third measurement was made at the end of 
the reach near the last sample site. Measurements were also made on 


major tributaries at the appropriate times. 


Discharges at each of the sample sites had diurnal fluctuations due 
to the operation of a hydroelectric plant upstream of Edmonton. Water 
levels at most sample sites were monitored every two hours so that the 
fluctuations could be matched with those from Water Survey of Canada 
water level gauges at Edmonton and Deer Creek in Saskatchewan. 
Fluctuations in discharge were estimated from winter stage-discharge 
curves for these sites. Discharges used in the following analysis were 
the estimated average discharges for the periods of interest. Table 2 
Summarizes the adopted discharges used in the travel time analysis. 
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Table 2 Summary of adopted reach-averaged discharges. 


Reach 


Edmonton Injection 
Edmonton - Ft. Saskatchewan 
Ft. Saskatchewan - Vinca 
Vinca - Waskatenau 
Waskatenau - Pakan 
Pakan - Shandro 


Shandro - Duvernay 


Shandro Injection 
Shandro - Duvernay 
Duvernay - Myrnam 
Myrnam - Elk Point 
Elk Point - Heinsburg 
Heinsburg - Lea Park 
Lea Park - Meridian 


Date 
Mar. 15 
Mar. 16 
Mar. 16 
Mar. 17 
Mar. 17 
Mar. 18 
Mar. 06 
Mar. 07 
Mar. 08 
Mar. 09 
Mar. 09 
Mar. 10 


16 


AV) 


18 


19 


07 
08 
09 


10 
1] 


10. 


Discharge 
(m/s) 


129 
108 
US 
145 
132 
153 


120 
124 
128 
134 
138 
142 


3.0 TRAVEL TIME ANALYSIS 
3.1 Channel Hydraulics 


The North Saskatchewan River between Edmonton and Meridian bridge 
has an irregular meander pattern with point bars, side bars, mid-channe] 
bars and occasional islands. Hydroelectric power production upstream of 
Edmonton produces a diurnal fluctuation of about 30% around the mean 
discharge during the winter. Winter flows also tend to be high because 
the river is regulated. The river is ice-covered from the middle of 
November to the middle of April except for a reach between Edmonton and 
Fort Saskatchewan which is ice free due to thermal effluent discharges 


from Edmonton. 


Ice cover formation on _ the North Saskatchewan River occurs by a 
combination of the juxtaposition of ice flows (which includes the 
redistribution of frazil) and by the consolidation of a juxtaposed cover 
in locations where the longitudinal forces on an advancing ice cover 
exceed the strength due to downward freezing of the juxtaposed cover. 
In the reach between Shandro and Duvernay, this consolidation process 


resulted in frazil accumulations approaching 2 m in thickness. 


Hydraulic characteristics were determined from cross-sections 
surveyed every five to ten kilometers over the whole reach. These 
Surveys were performed over the month of January by crews from Hamilton 
and Olsen Surveys Ltd. The ice thickness was found to vary between 
0.3 m - 0.8 m while the river top width varied between 120 m and 330 m. 


The average times of travel of the leading edge, the peak, the 
centroid, and the trailing edge of the dye cloud for each subreach are 
listed in table 3. As well, the cumulative travel times for each test 
are plotted against distance from injection in figures 2a-b. The slope 
of the lines in these figures are inversely proportional to the 
velocities. For example, in both cases the velocity of the peak is 
about 5% faster than the velocity of the centroid. 
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Table 3 Subreach average travel times. 


Subreach Length Travel Times 
Leading Peak Centroid Trailing 
Edge Edge 
(km) (hrs) (hrs) (hrs) (hrs) 
Edmonton 
-Fort Saskatchewan 36.2 14.48 15.50 16.75 21.86 
Fort Saskatchewan 
-Vinca AS) 5) 11.78 14.50 16.24 29.24 
Vinca 
-Waskatenau 25.6 SROs T3038 14.17 ies5 
Waskatenau 
-Pakan 28.1 Ze 4a We) SAS) 7/ 14.97 25.62 
Pakan 
-Shandro (Alp a2 12.50 14.02 14.30 13.56 
Shandro 
-Duvernay ilaer3 26.22 28.98 29-2 30.46 
Shandro 
-Duvernay 51.4 26.96 29.00 31.56 40.37 
Duvernay 
-Myrnam SS 14.60 18.03 18.18 (60 (42 
Myrnam 
-Elk Point 25538 16.27 14.47 16.02 19.74 
Elk Point 
-Heinsburg 33.4 15.89 18.53 18.82 (fae VY) 
Heinsburg 
-Lea Park 18.4 102.10 13.97 12.82 15.53 
Lea Park 


-Meridian 21) <S 15 92 14.00 15293 20.08 
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Figure 2a Growth of travel times with distance from injection, Edmonton 
to Duvernay. 
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Figure 2b Growth of travel times with distance from injection, Shandro 
to Meridian. 
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This difference in velocity between the peak and the centroid of 
the time-concentration curves raises the question of which quantity 
represents the average velocity of the water. Obviously, the centroid 
is the measure of the average velocity of the dye mass but this velocity 
only corresponds to the average water velocity in a uniform laboratory 
channel where all areas of the channel contribute to the flow. In 
natural streams, however, especially for ice-covered winter conditions, 
the cross-section variability causes zones of stagnant or very slow 
moving water. These zones contribute little if anything to the 
discharge therefore they should not be included in the determination of 
average velocity. The travel time of the centroid of the dye cloud is 
affected by these stagnant zones especially if they produce long 
drawn-out receding limbs. Alternat ively, the cross-section average 
velocity of the peak concentration incorporates the cross-sectional 
variations in velocity but does not include the effects of the storage 
and release of dye from stagnant areas. Therefore, the velocity of the 
peak should be used to calculate the hydraulic characteristics of a 


river reach. 


Generalization of the hydraulic characteristics measured at the 
individual sites into reach-averaged characteristics for each of the 
Subreaches was done by defining the mean velocity, U for each subreach 
as x/t where x 1S the reach length and an is the time it takes for the 
peak of the dye cloud to pass through the reach. The resultant subreach 


mean velocities are given in table 4. 


Using the time of travel of the peak to determine the mean velocity 
may lead to underestimating the time of travel if mixing across the 
channels not fully established. That is, the dye moves primarily in the 
zone of high velocity and, hence, the velocity is over estimated because 
the entire flow is not represented. To estimate the magnitude of this 
measurement error in the initial subreach, the time of travel between 
Shandro and Duvernay was measured twice, once at the beginning of a test 
reach and once at the end of a test reach. Times of travel could not be 
compared directly because the discharges were not the same for both 


tests; however, velocities could be compared assuming the channel 


M5 . 


Table 4 Summary of average hydraulic characteristics in late winter 


for each subreach. 


Subreach Reach Vel-  ODis- Slope Top Flow Composite 
Length ocity charge Width Depth Roughness 
(km) (m/s) (m/s) (m)  (m) 
Edmonton 


-Ft. Saskatchewan 365 5 Ob65 Us) DO0CS66 12) sale 0.032. 


Ft. Saskatchewan 


-Vinca 29.57; 0.57 108 0.000381 198 0.99 0.022 
Vinca 

-Waskatenau 25.6 0.56 115 0.000443 232 0.90 0.023 
Waskatenau 

-Pakan 28.1 0.56 145 0.000382 145 1.28 0.026 
Pakan 

-Shandro 212 0.54 132 0.000397 187 1.30 0.028 
Shandro 

-Duvernay 51.3 0.49 153 0.000359 218 1.43 0.031 
Shandro 

-Duvernay 51.4 0.49 120 0.000359 218 jegoaly2 0.026 
Duvernay 

-Myrnam SS 42 0.51 124 0.000216 210 ue as) 0.020 
Myrnam 

-Elk Point 25.3 0.49 128 0.000220 264 1.00 0.019 
Elk Point 

-Heinsburg 33.4 0.50 134 0.000464 212 Leo Zy 0.032 
He insburg 

-Lea Park 18.4 0.37 138 0.000232 214 1.74 0.038 
Lea Park 

-Meridian , 27.9 0.55 142 0.000179 246 1.01 0.016 


* 
open water reach 
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roughness remained relatively constant. The velocity in the initial 
Subreach was found to be overestimated by about 9%. The velocity in the 
initial subreach between Edmonton and Fort Saskatchewan is likely 
overestimated by less than this value because three, rather than two, 
injection points were used thereby spreading the dye more quickly over 
the channel width. 


The mean velocity and the adopted discharge, Q were used to 
determine the flow area, A, for each Subreach. From the flow area and 
top width, the mean depth, d, and the hydraulic radius, A/2W, could be 
calculated. From these, the composite Manning roughness, no» was 
determined by: 

[1] CCAM, AS. 


° 


)/Q 


where S is the average water surface slope. Table 4 summarizes the 
results. 


The composite Manning roughness was found to vary from 0.015 to 
0.038. This range of roughness is similar to the open water roughnesses 
given by Kellerhals, Neill and Bray (1972). This suggests that the 
underside of the ice cover is very smooth and does not affect the 
composite roughness significantly. When a similar analysis was done on 
the Athabasca River between Hinton and Athabasca (Andres, Van Der Vinne, 
and Trevor, 1989), no significant increase in roughness was observed for 
low flow winter conditions; however, studies done on the Wapiti and 
Smoky rivers indicated a significant increase in roughness due to the 
ice cover in some areas. This difference in the effects of an ice cover 
on the hydraulic roughness cannot be explained at this time. More 
effort must be made to quantify the ice cover characteristics if this 
difference is to be resolved. 


3.2 Generalized Travel Times 


The travel times for any given discharge are related to the mean 


velocity of the flow. Because the mean velocity varies monotonically 
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with discharge in most cases, the measured travel times or velocities 
are only valid for the discharges at which the measurements were carried 


out. 


The extrapolation to other discharges can be made by combining the 


equation of continuity: 
[2] U = Q/Wd 


with a roughness equation such as the Manning equation: 


[3] Q = (1/n,) R Wd S 


for which all the variables have been defined previously. 


It should be noted that for a wide, ice covered channel, the 
hydraulic radius, R, 1S given by: 


[4] R = A/2W = d/2 


to account for both the top and bottom boundaries. Combining equations 
[2], [3], and [4] produces equation [5] which defines the mean velocity 
in terms of the width, the slope and the composite roughness 


[5] Wee(O/aw) «(So fn) 


Thus for a given river reach, with known values of S, W and n,» the 
velocity, and hence the travel time, can be obtained for any given 
discharge. Figure 3 shows equation [5] in graphical form along with the 
test results. The curves show the general trend of increasing velocity 
with increasing discharge for various values of slope and roughness. 
The data fits equation [5] exactly because a form of the equation was 
used to calculate the composite roughness. 


Tor: 


Ice cover 
Open water 


VS/n =1.5 


V3/n5=1.0 


VS/n_=0.5 
(°) 


Velocity (m/s) 


VS/n =0.25 


Discharge per unit perimeter (m/s) 


Figure 3 Relationships between mean velocity and discharge per unit 
width for the North Saskatchewan River under an ice cover. 


The above analysis assumes that the composite roughness is constant 
for all flows, however, it can increase significantly with decreasing 
discharge as the flow depth decreases in relation to the variations in 
the channel bottom. If this relationship is approximated by 


[6] noe a). 


then equation [5] becomes 


[7] U=C Oe eae”. 


where C is a function of the width and slope of the channel. In some 
cases, it may be necessary to perform additional travel time 
measurements at other discharges to define the exponent b. On the North 
Saskatchewan River, variations in roughness due to changes in discharge 


19. 


should be minor because the river is_ regulated. However, large 
variations in roughness can occur from year to year, depending on the 
hydrologic and meteorologic conditions. The impact that the ice cover 
has on the composite roughness is not well understood and therefore it 
1S recommended that the data not be extrapolated too far beyond the 


conditions for which the measurements were made. 
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Table 5 Reach-average hydraulic characteristics. 


Reach Discharge Slope Width Velocity Depth Shear Composite 
3 Velocity Roughness 
(m /s) (m) (m/s) (m) (m/s) 
Edmonton- 
Duvernay 35 0.000389 196 O57 WZ 05048 0025 
Shandro- 
-Mer idian 135 0.000294 225 0.49 he22- 02024 0.025 


4.0 DISPERSION ANALYSIS 


The tracer dye studies which were done to measure the hydraulic 
characteristics of the North Saskatchewan River can also be used to 
measure the rate of longitudinal dispersion, the process by which 
concentrations are reduced as the pollutant cloud spreads along the 
length of the channel. This longitudinal dispersion rate can be used to 
model or predict the extent of toxic concentrations downstream from an 


accidental chemical spill. 
4.1 Reach-average Hydraulic Characteristics 


Reach-average hydraulic characteristics are required to determine 
the nature of longitudinal dispersion in the North Saskatchewan River 
from the data obtained in the tracer dye tests. Reach-average 
velocities can be obtained from the slope of the time-to-peak versus 
distance data for each reach as seen in figure 2. Velocities were quite 
consistent in both test reaches. Reach-average values of slope and 
width were calculated from the subreach values given in table 3. 
Reach-average flow depths were calculated from the average discharge, 
velocity and widths for each reach. All these hydraulic characteristics 
are given in table 5. 
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Table 6 Summary of mixing lengths for single centerline injection. 


Reach Reach Vertical Transverse Est imated 
Length Mixing Mixing Linear Mixing 
Length Length Length, a, =1.0 
(km) (m) (km) (km) 
Edmonton- 
Duvernay 198.1 103.6 142 378 
Shandro- 
Meridian 189.5 (A5) 6 Tf 182 486 


4.2 Mixing Zones 
4.2.1 Vertical Mixing 


The mixing process in rivers begins as a three-dimensional process 
but the pollutant concentration quickly becomes uniform over the depth. 
According to Elhadi et al (1984), vertical variations in concentration 
are less than 2% at a distance L, downstream from a source located at 
mid-depth when 


Co 18 d U/U 
where U. is the shear velocity defined by 
1/2 

[9] U, = (QRS) 
in which g 1s the acceleration due to gravity. 

At Shandro, where the injection was at mid depth, the vertical 
mixing length was about 26 m (table 6). At Edmonton, however, the 
injection was near the surface so the mixing length should be four times 


greater than that given by equation [8], that is, 104 m. Thus, it can 
be assumed that vertical mixing is instantaneous compared to the total 


ie ee i ee ee el eee ee ed = 


LL. 
reach lengths. 
4.2.2 Transverse Mixing 


Once the pollutant concentration becomes uniform over the depth, 
the mixing can be treated as a two-dimensional problem. Further 
downstream, the pollutant concentration becomes essentially uniform over 
the width of the channel as well. Transverse mixing is essentially 
complete at a distance La from a source located at mid-width when 


[10] L, = (0.1U/e,)(W) 


(Fischer et al, 1979) where e, is the transverse mixing coefficient 


usually estimated from 
aya e, = k,dUu, 
When equations 10 and 11] are combined the result is 


[12] L, = (O.1/k,) (U/U,) (W°/d) 
Yotsukura and Cobb (1972) introduced a parameter, Ko which varies with 
the position of the source in the section and the degree of mixing 


[13] Le LG Roy UR cei 


Values of this position parameter for various fractions of cumulative 
discharge across a section can be obtained from table 7. 


Values of k, range between 0.4 and 0.8 for natural streams (Fischer 
et al, 1979); however, these values were measured in open water. 
Engmann (1974) found that in laboratory channels values of e_ for flows 
with top covers were only half as large as those for corresponding free 
surface flows. Engmann proposed that €, was proportional to the 
hydraulic radius, R rather than the flow depth, d. Beltaos (1978a), 


however, measured transverse mixing under ice covers in natural streams 
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Table 7 Position parameters for various fractions of cumulative 
discharge across a section and 95% complete mixing. 


Fraction of Cumulative Position 
Discharge from Nearest Bank Parameter 


(on i an i ep il ap i an i ae) 
ot & G Por © 


and found that e, was more a function of d than R. This conflict has 
not been resolved so it best to use the field data as a guide rather 
than the laboratory data. Beltaos (1979) also measured e, in the North 


Saskatchewan River without an ice cover and found k, to be 0.17. 


The transverse mixing lengths listed in table 5 were calculated 
using equation 13 with KS = 15.7 which corresponds to a single injection 
in the center of the channel. These lengths define the points where 
mixing is 95% complete. In the actual tests, however, multiple 
injections were used to reduce the length required for transverse 


mixing. 


At Shandro, two injections were used, each at one-third of the 
channel width. This corresponds to 0.2 and 0.6 of the cumulative 
discharge distribution from the left bank and, from table 7, values for 
kK, of 4.5 and 7.2 respectively. Transverse mixing lengths for these 
injections individually are 635 km and 397 km respectively, well beyond 
the 189.6 km length of the test reach. Complete mixing to the nearest 
bank from each injection, however, will occur over a much shorter 
distance. Still, the left bank concentration will be higher than the 
right bank because the discharge in the left zone was less than that for 
the right zone. After this point, the only transverse mixing that wil] 


occur will be the leveling of the concentrations between the two zones. 


To calculate these shorter mixing lengths, it is assumed that the 
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Table 8 Transverse mixing lengths for multiple injections. 


Location Cumulative Position Total Effective Mixing 
Discharge Parameter Transverse Width Length to 
from nearest Mixing Length nearest bank 

bank (km) (m) (km) 

EDMONTON 

Left side 0.4 Tine 308 160 94 

Centerline 0.3 S53 384 120 53 

Right side Oot Bg 567 40 Joe 

SHANDRO 

Left side OFZ 4.5 635 90 29 

Right side 0.4 Te2 396 180 116 


injections occur at the center of a streamtube with a width twice the 
distance from injection point to the nearest bank. For q/Q = 0.2, this 
WiidiEMe VS cl 0.2225) =, 90m and, for, g/Q =,0.6,. the width 1s. 2(0.4*225) = 
180 m. Using these widths and kK, = 15.7 for a centerline injection, the 
transverse mixing lengths become 29 km and 116 km respectively 
(table 8). Hence, after 116 km, the river has two completely mixed 
concentration zones which are slowly equalizing. This can be seen in 
the measured concentration distributions given in Appendix A. At Elk 
Point, 110 km downstream of the injection, the right side of the channel 
still has a significantly lower concentration than the left side. 
However , at Heinsburg, 143 km downstream of the injection, the 


concentrations across the channel are much more uniform. 


At Edmonton, three injection points dividing the width uniformly 
were used to reduce the initial concentrations. The discharge 
distribution at the Edmonton injection site is unknown because the 
discharge was measured downstream at the Capilano bridge rather than at 
the injection site. However, the distribution can be estimated assuming 
that the cross-section is triangular with the deepest point near the 
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left bank on the outside of the bend and that the local average velocity 
is directly proportional to the depth. Also, it was observed during the 
injection that the inside or shallow third of the width was blocked by 


Vee: These assumption result in cumulative discharges of 0.4, 0.7, and 
OFS trom the, Vert” bank. The transverse mixing lengths for these 
injections are given in table 8. Total transverse mixing is not 


completed in the 198.1 km test reach; however, within 100 km the dye is 
fully mixed to the banks. There is still a small] concentration gradient 
from one bank to the other after this distance but it can be ignored if 


cross-section average values are used in the analysis. 
4.2.3 Linear Dispersion 


When transverse mixing is complete, further mixing is due to 
longitudinal dispersion, the process which increases the spread of 
pollutant along the channel length as the pollutant travels downstream. 
Data from a large number of field studies by numerous investigators 
indicate that the onset of Fickian dispersion is delayed for some 
distance downstream after transverse mixing is completed (Beltaos, 
1978b). In this linear dispersion zone, differential advection and 
variations in channel geometry cause longitudinal dispersion to occur 
more quickly than predicted by the Fickian dispersion model discussed in 
the following section. The field measurements suggest that in this zone 


the standard deviation, o, of the time-concentration distribution grows 


t 
linearly with distance; hence, this type of dispersion is called 


‘linear’ dispersion. 


The location of the transition from linear dispersion to Fickian 


dispersion can be estimated using 
[14] L, = a, (U/U,)(W/d) 


where a, is a factor varying between 0.48 and 1.8 depending on the 


L 
degree of channel irregularity. The transition occurs between L, and 


3L,. This linear mixing length can be estimated for single injections 


by assuming that a, = 1.0. Estimated values of L, are given in table 6 
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along with the other single injection mixing lengths. 


Approximate values of linear mixing lengths for multiple injections 
can be obtained by assuming each injection stays within its own 
streamtube. Hence, for two injection points the effective width is 
one-half of the river width and for three injection points the effective 
width is one-third of the river width. These multiple injection linear 
mixing lengths are much shorter than those for single injections because 
the mixing lengths are proportional to width squared. Using this 
procedure, the estimated value of L, for the Edmonton injection is 42 km 


while for the Shandro injection it is 121 km. 


No analytical solution exists for linear dispersion so Beltaos 
(1978b) proposed the following empirical solution based on the existing 
data 


1 
[15] Geel [= exp|1 -= I) 
where 6 1S a parameter defined by 


Or 

le] 6 = — 
U/x 

The calculation of 6 from measured concentrations can be simplified 

by defining the pollutant spread in terms of the half-duration, AT 


rather than the standard deviation, o This eliminates the problems in 


t: 


evaluating o, due to incomplete measurements on the falling limb. The 


t 
half-duration is defined as the period of time during which the 
concentration is greater than one-half of the peak concentration 
(figure 4). The half-duration, AT of the empirical curve described by 
equation [15] is equal to 2.36 o,; therefore, B can be defined in terms 


of AT 


[17] Bie 0.13 [At] 
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Concentration 


Time from injection 


Figure 4 Definition sketch for variables describing time-concentration 


distributions 


4.2.4 Fickian Dispersion 


The final mixing process is termed Fickian dispersion. This 


process can be modeled using the one-dimensional mixing equation 


2 
[18] 2 ele falls 
BE won| bn 


which, for a mass, M _- instantaneously released throughout a 


cross-section, has the solution 


2 
[19] pe eee exp| ees) 
2QVaDt 4Dt 


where C is the cross-sectional average concentration. The properties of 
this analytical solution indicate that the dispersion coefficient, D can 


CS). 


be defined by 


[20] D = — 


20 : 
That is, the rate of change of the variance, o, with distance is 


constant. 
4.3 Dispersion Parameters 


The dispersion parameters discussed in the previous sections can be 
evaluated using variables describing the -cross-sectional average 
time-concentration distributions. The variables describing these 
Ole at each location are: time-to-peak, th half-duration, AT; 
variance, o, ; peak concentration, Ch: and mass recovery ratio, M/M.. A 


summary of these variables for each location is given in table 9. 


The linear dispersion parameter, 6 can be evaluated using equation 
[17], that is, reach-average values of £6 for the North Saskatchewan 
River can be obtained from the slopes of the best-fit lines through the 
half-duration versus time-to-peak data shown in figure 5. Values of £B 
for the two test reaches are listed in table 10. These values of £6 are 
lower than values measured in other ice covered rivers; however, they 
are consistent with the trend of increasing 6 with U./U shown in 
figure 6 (Beltaos, 1978b). The data shown in this figure also suggests 
that the presence of an ice cover tends to increase the magnitude of the 


linear dispersion parameter. 


Values of D can also be evaluated from dye test data provided that 
the test reaches extend into the Fickian range. As stated in equation 
L20]}. hs 
the slope of this line. Three data points are needed to establish that 


vs x iS linear in the Fickian range and D is proportional to 


the trend is linear. The values of D calculated from slopes obtained 
from figure 7 are given in table 10. There are some significant 
deviations from the best-fit line for the Edmonton injection, likely due 


to the lengths of the recession limbs included in the variance 
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Table 9 


distribution variables. 


Location 


Edmonton injection 
Ft. Saskatchewan 
Vinca 
Waskatenau 
Pakan 
Shandro 
Duvernay 


Shandro injection 
Duvernay 
Myrnam 
Elk Point 
Heinsburg 
Lea Park 
Mer idian 


Distance 
from 
Injection 
(km) 


S19) 
Gor 
Oil 
119. 
146. 
1987 


— 0O D 1 WH MO 


ile: 
84. 
109: 
143. 
161". 
che) 
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Table 10 Summary of dispersion parameters. 


Reach 


Edmonton- 


Duvernay 


Shandro- 


Meridian 


Linear 
Dispersion 
Parameter 


B 


0.00612 


0.00429 


Time Half - Variance 
to duration 
Peak 
(hrs) (hrs) — (hrs~) 
155 LoS 5 Ad 
S00 6.20 9.56 
43.0 9.26 10e23 
57 50 1134 29.94 
Ne 18} 69) 32.47 
100.0 2 45.83 
29.0 4.95 4.61 
47.0 8.08 Hel 
Oia 9.39 16.33 
80.0 13.16 27213 
94.0 14.98 37.75 
108.0 15.24 47.38 
Linear Linear 
Mixing Mixing Length 
Length Factor 
Ly Me 
(km) 
34.6 0.83 
66.4 O55 


P 


( 


OoOoOrr fh oO 


me AIP W LP 


Cor 


Summary of cross-section average time-concentration 


eak Mass 


Concen- Recovery 


tration Ratio 


ug/L) 
.875 0.750 
ON 0.748 
.298 0.756 
.127 0.849 
.902 0.824 
.676 0.832 
.973 0.783 
.584 0.934 
. 843 0.899 
.085 0.922 
811 0.953 
. 753 0.966 
Fickian 
Dispersion 
Coefficient 
D 
(m/s) 
344 


(E5)3I 


ge rs 


ee i ae.) | 
yo! ipaangae 
 .aetaatie 
‘“y208 | ~ Meet sanctwe¥ ish ‘ontT “waged 
yvaoeN ee nottewe of 


 @ 
A A 


Ff 
Bit) Nt 


7 my mont ; | : | , Were mi : ie ry eat 
ot tet nok iar! | 7 ae 4504 nottoetnl | ~ 


| Pe ae. lee 
(toy (em) a oe? ey ee 


a. 


Saat Re. 
os 
~~ 
© 


027.6 ave.d OS.f ett 
Oot.0 ©=6— $8.5. Wee 0S.a 
“oay 0 eeg.t €5.0F  @6.e 
Gas.0 {S1.f. © 80.85 $f .04 
$8.0 9oe.0 8 Vh.gt ef .ef 
$e2.0 av9.0° E89. dp Sf.Mf 


ote 47 


agi azepat aban 


cay .o cve.- tae 22.5 0,e5 bile veo 

bf? .o eBs.& «th .' 80.8 0,14 a, 58 | I ee 
eee. co6.& FE.9! €t.@ aie ee. Intee BF 
gse.0 o80.5 ef.¥3 a[.=! 5.0e €. chi i pudanigh 
cée.0 (18.5 e.ve Be.ti 0.6e t. fat | at eet 7 

aee.0 86a. Bf. &$.25 0,902 ‘dé, 08: | teothiaeh 


. 2 19S ome T8q poreroqate Yo yrsomez Of sidet | 
aetag snonts saont J wets "| damea 7 . 
aot zreqett | d3peeas gntzih = garci® not ereqe td a eee ¥ a 
ine! 91 Yeo) voles? (‘ie otometed 
ee i ; i? © < > Le 1 es 
7 “0 f fe a y" a 7 4 im, : =. oe 
{e\*n) ae i a ; 


a 


s 


30. 


Edmonton Injection 
Shandro Injection 
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Figure 5 Variation of half-duration with time-to-peak from injection. 
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Figure 6 Variation of 8 for both open water and ice covered conditions 
(Beltaos, 1978b). 
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calculations. The dispersion coefficients measured in this study are 
consistent with values reported in the literature for channels of this 
width and roughness (figure 8). Not enough of these measurements have 
been done under an ice cover to make any conclusions regarding the 


effects of an ice cover on the dispersion coefficient. 


Actual values, or at least lower limits, of L, can also be obtained 
from figure /7. Beltaos (1978b) defined L, as the x-intercept of a 
straight line fitted through the o,- 
zone. Once L, is evaluated, the actual value of a can be obtained from 


vs x data in the Fickian dispersion 


equation [14]. Values of both L, and a, for the test reaches are given 
in table ~10. The three injection points at Edmonton cause Fickian 
dispersion to begin much more quickly in the upper test reach. In the 
lower reach, the reach length was slightly less than 3L, away from the 


injection site so Fickian dispersion did not occur. 


The values of L, given in table 10 are similar to the values 
estimated in section 4.2.3 but are less than the transverse mixing 
lengths given in table 8 because L, is evaluated from cross-section 
average values. Thus, only average concentrations can be predicted 


before transverse mixing is complete. 
4.4 Dispersion Model Predictions 
4.4.1 Linear Dispersion 


Time-concentration distributions were generated for each sample 
Site using the linear dispersion model defined in equation [15]. These 
predicted distributions are plotted in figure 9 along with the 
cross-section average of the measured distributions. The major 
discrepancy between the measurements and the model predictions is in the 
magnitude oof the peak concentrations. The model predicts peak 


concentrations as much as 50% higher than the measured values. 


An assumption in the model is that the dye is completely conserved 
but the mass recovery ratios plotted in figure 10 show that as much as 
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Figure 7 Increase in variance with distance from injection. 
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Figure 8 Variation of dimensionless dispersion coefficient 
(Beltaos, 1978b). 
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Figure 9 Comparison of measured time-concentration distributions with 
those predicted by the linear dispersion model. 
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Mass recovery ratio, M/M 
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Figure 10 Variation in mass recovery ratio with distance from injection. 


25% of the dye was undetected at some of the sample sites. This 
undetected dye may be lost due to chemical and physical processes or it 
may be due to inaccurate assumptions of discharge or water temperature. 
For example, the increase in mass recovery ratio between Fort 
Saskatchewan and Pakan is likely because water temperatures in the open 
water reach were greater than the assumed near freezing value. If the 
model were modified to account for losses in mass, the predicted 
time-concentration distributions would match the measurements much more 
closely. 


The peak concentrations predicted by the model may also be higher 
than the measured peaks because of the difference in shape between of 
the predicted and measured time-concentration distributions shown in 
figure 9. The more symmetrical curves predicted by the model have 
shorter durations, therefore the peaks would be higher even if the dye 
masses were the same. The model could be improved by making it predict 
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the shape of the time-concentration distributions more accurately. 


The times-to-peak of the predicted distributions are slightly 
different from the measured values as well. This difference is caused 
by the model using reach-average velocities rather than measured 


Sub-reach velocities. 


The model predictions are reasonable considering all the 
assumptions and approximations used to obtain the modeling parameters. 
The linear dispersion model can be _ used with confidence to predict 
concentrations downstream of accidental releases of pollutants as long 
as some margin for error is incorporated into the results. However, 
this model should only be used until the point where Fickian dispersion 


begins. 
4.4.2 Fickian Dispersion 


Elhadi et al (1984) state that once a pollutant cloud passes a 
point 3L, away from the spill site, equation [19] should be used to 
predict concentrations. Equation [14] can be used to estimate be for 
Single injections; however, this distance is so large for the North 
Saskatchewan River that Fickian dispersion would not occur until well 


into Saskatchewan. 


In the case of the multiple injection tests done in this study, L, 
can be estimated by dividing the channel width by the number of 
injections and using this reduced width in equation [14]. As well, 
Equation [19] cannot be used as it stands because its derivation assumed 
complete mixing across the section immediately upon injection. To make 
equation [19] suitable for modeling dispersion from sources other than 
completely mixed sources, it must be assumed that the pollutant is 
injected, fully mixed, at a virtual origin a distance L, downstream from 


the actual injection. Thus equation [19] becomes 
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where be 18 defined as Bue 


Predicted distributions using this modified Fickian dispersion 
model] are shown in figure 9 along with those of the linear model. All 
locations’ Gpeatem than a were modeled to show the inaccuracy of the 
model in the linear dispersion zone. Note that the peaks predicted by 
the Fickian dispersion model are highest at both Vinca and Myrnam; 
however, by the end of the reaches, predictions from the two model are 
quite similar. If the reaches were extended, the linear model would 
begin to overpredict dispersion while the Fickian model would be more 


accurate. 


Se 
4.0 SUMMARY 


Two tracer dye tests were carried out on the North Saskatchewan 
River between Edmonton and the Saskatchewan border. Discharges during 
the tests ranged from 108 m/s OpLoS m/s. Ice thicknesses varied 
between 0.3 m and 0.8 m with frazil accumulations approaching 2 m in 
some places. Velocities calculated from the travel times of the peak 
concentration ranged between 0.37 m/s and 0.57 m/s in the ice-covered 
reaches. Composite hydraulic roughnesses varied between 0.016 and 
0.038. The extrapolation of these velocities is possible, provided some 
confidence can be placed on the adopted values of hydraulic roughness 
for the extrapolated conditions. Because the roughness depends 
Significantly on, the ice conditions, which in turn can vary 
substantially from year to year and throughout the winter season, it is 
recommended that care be taken to verify the characteristics of the ice 


cover prior to extensive extrapolations. 


Two types of longitudinal dispersion occurred during the dye tests. 
First, after transverse mixing was completed, ‘linear’ dispersion 
occurred, where the length of the dye cloud grew linearly with distance. 
After a transition period, the second type of longitudinal dispersion, 
Fickian dispersion, began where the length grew with the square root of 


distance. 


The linear dispersion parameter was found to be 0.00612 for the 
Edmonton injection and 0.0429 for the Shandro injection. These values 
are lower than the values measured on other ice covered rivers but this 
is due to the low values of hydraulic roughness found in these river 
reaches. The Fickian dispersion coefficient was found to be 344 m/s 
for the upper reach and 291 m/s for the lower reach. These values are 
consistent with values reported in the literature for channels of 


Similar width and roughness even with the presence of an ice cover. 


Once the loss of dye was taken into account, the major discrepancy 
between the results of the longitudinal dispersion models and the 
measurements was the shape of the concentration distributions. The 


models predicted much more symmetrical distributions than were actually 


38. 


observed. 


The analysis of the data presented herein provided some information 
on both the hydraulic and mixing characteristics on the North 
Saskatchewan River; but, further information 1s required to gain a more 
complete understanding of these characteristics. llomigather this 


information, it 1S recommended that: 


(1) a low flow open water travel time test be undertaken to better 
quantify the low flow open water hydraulic and mixing 


characteristics and 


(2) a study be undertaken to quantify the variations in length of 
the open water reach downstream of Edmonton over winter to 
determine the impact of this variation on the hydraulic and 


mixing characteristics. 
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Dye Concentration Measurements 
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Concentration (ug/L) 


Depth (m) 


Duvernay sample site 
2.315km downstream of Hwy 36 bridge 
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Center Line 
Right Bank 
Average 
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Concentration (pg/L) 
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Myrnam sample site 
at old ferry access 620m downstream of Hwy 881 brigde 
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Concentration (ug/L) 
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Elk Point sample site 
125m upstream of Hwy 41 bridge 
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Center Line 
Right Bank 
Average 
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Concentration (ug/L) 


Depth (m) 


Heinsburg sample site 
750m downstream of Hwy 893 bridge 
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Center Line 
Right Bank 
Average 
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Concentration (ug/L) 


Depth (m) 


Lea Park sample site 
30m upstream of Hwy 897 bridge 


Left Bank 
Center Line 
Right Bank 
Average 
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Concentration (g/L) 


Depth (m) 
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Meridian sample site 
30m downstream of Hwy 17 bridge 
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Left Bank 
Center Line 
Right Bank 
Average 
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